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ABSTRACT 

I present high-resolution column density maps of two molecular clouds having strikingly different 
star formation rates. To better understand the unusual, massive G216-2.5, a molecular cloud with 
no massive star formation, the distribution of its molecular gas is compared to that of the Rosette 
Molecular Cloud. Far-infrared data from Herschel are used to derive A(H 2 ) maps of each cloud and 
are combined with ico data to determine the CO-to-H 2 ratio, Aco- In addition, the probability 
distribution functions (PDFs) and cumulative mass fractions of the clouds are compared. For G216- 
2.5, (iV(H 2 )) = 7.8 x 10 20 cm” 2 and (Xco) = 2.2 x 10 20 cm” 2 (K km s” 1 )” 1 ; for the Rosette, 
(A(H 2 )) = 1.8 x 10 21 cm” 2 and (Aco) = 2.8 x 10 20 cm” 2 (K km s” 1 )” 1 . The PDFs of both clouds 
are log-normal for extinctions below ~ 2 mag and both show departures from log-normality at high 
extinctions. Although it is the less-massive cloud, the Rosette has a higher fraction of its mass in 
the form of dense gas and contains 1389M 0 of gas above the so-called extinction threshold for star 
formation, Ay = 7.3 mag. The G216-2.5 cloud has 874M 0 of dense gas above this threshold. 

Subject headings: ISM: clouds — ISM: individual objects (G216-2.5) — ISM: individual objects 
(Rosette) — ISM: dust, extinction infrared: ISM — submillimeter: ISM 


1. INTRODUCTION 

Since stars are observed to form in molecular clouds 
(MCs), the properties of these stellar nurseries are ex¬ 
pected to be closely tied to the processes of star forma¬ 
tion. It has been a continuing goal of astrophysicists 
to thoroughly characterize the structure of MCs and to 
understand how this structure relates to protostellar col¬ 
lapse. The realization since at least the 1980s that MCs 
of similar masses can have drastically different rates of 
star formation gives further impetus to fulfilling this ob¬ 
jective. 

The G216-2.5 Molecular Cloud and the Rosette Molec¬ 
ular Cloud, as emblems of extreme contrasts in star 
formation activity, are good astrophysical laboratories 
for testing issues related to MC and star formation. 
Also known as Maddalena’s Cloud or the Maddalena- 
Thaddeus Cloud, G216-2.5 is noteworthy for its low tem¬ 
perature (~ 10 K), high mass (~ 1O 5 M 0 ), but low lev¬ 
els of star formation activity, compared to other MCs 
(Maddalena & Thaddeus 1985). In a series of papers 
in the 1990s, Lee, Snell & Dickman (1991; 1994; 1996) 
mapped G216-2.5 in 12 CO and 13 CO and studied its in¬ 
frared properties. Lee et al. confirmed the earlier re¬ 
sult of Maddalena & Thaddeus (1985) that the molecu¬ 
lar cloud does not host massive star formation. Through 
near-infrared imaging of point sources detected by the In¬ 
frared Astronomy Satellite (IRAS), they suggested that 
low-mass star formation is occurring primarily toward 
the boundaries of the cloud, and deduced that it is a rem¬ 
nant cloud that achieved its properties from past episodes 
of massive star formation. Only recently were observa¬ 
tions presented that indicate low-mass star formation is 
occurring in the main body of the cloud (Megeatlr et 
al. 2009). The Rosette Molecular Cloud stands in stark 
contrast to G216-2.5. The well-studied MC contains the 
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OB association NGC 2244 whose stars excite the Rosette 
Nebula which, in turn, may be triggering star formation 
elsewhere in the cloud (Heyer et al. 2006). Mid-infrared, 
near-infrared, and X-ray observations have revealed the 
presence of several embedded clusters within the Rosette 
(e.g., Phelps & Lada 1997; Roman-Zuniga et al. 2008; 
Ybarra et al. 2013). 

Thus, the Rosette and G216-2.5 represent two limits 
of star formation activity: on the one hand, the star 
formation rate (SFR) per unit mass in the Rosette is 
higher than the Galactic average (Williams & McKee 
1997), while G216-2.5 does not contain even one O star, 
as is expected for most Galactic MCs of similar mass. 

Does G216-2.5 have such low levels of star formation 
because it is very young (a few Myr) and has only begun 
to form stars, as suggested by Maddalena & Thaddeus 
(1985)? Or is the MC in between episodes of massive 
star formation, as Lee et al. (1994) proposed? The main 
disadvantage of the latter proposal is the lack of stars 
detected in G216-2.5. A study of the atomic hydrogen 
(H I) gas associated with the cloud was more ambigu¬ 
ous in its conclusions about the cloud’s evolutionary sta¬ 
tus (Williams & Maddalena 1996). Moreover, the argu¬ 
ments that G216-2.5 is a relic of past star formation were 
made prior to the Spitzer observations by Megeath et al. 
(2009), which show there is an embedded population of 
young stellar objects (YSOs) in the main body of the 
cloud. 

In this and a following paper, I take a fresh look at 
G216-2.5 to address the questions of the MC’s evolu¬ 
tionary status and why it has such a low star formation 
activity. In light of new observations since the 1990s, es¬ 
pecially those of the Herschel Space Observatory, as well 
as advances in numerical simulations of MC formation 
and evolution, a reexamination of this unusual molecular 
cloud is timely. This paper addresses the G216-2.5 mys¬ 
tery by comparing properties of its molecular gas- in 
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particular, the column density distribution -with that of 
the well-studied Rosette Molecular Cloud, which is con¬ 
sidered to be representative of MCs harboring high-mass 
star formation. In §2, I review previous results that are 
germane to this study. In §3 and §4, I discuss the data 
and methods used in this study. The results, including 
column density maps and CO-to-H 2 conversion maps, are 
presented in §5. The key findings are further discussed 
in §6 and summarized §7. 

2. BACKGROUND 

In the 1990s, Lee et al. (1991; 1994; 1996) reported 
on a rigorous set of investigations into the molecular 
gas, star formation, and dust properties in G216-2.5. 
Using the FCRAO 14 nr telescope, Lee et al. (1991) 
mapped G216-2.5 in the 12 CO and 13 CO J = 1 — 0 tran¬ 
sitions. They showed that with 12 CO antenna tempera¬ 
tures of less than 4 K across the cloud, G216-2.5 is colder 
than most observed MCs. In the second paper, Lee et 
al. (1994) described how the large 12 CO line widths, 
~ 8 km s -1 , might have resulted from prior star forma¬ 
tion. Their data indicate a complex velocity structure 
(e.g., shells and rings) that could be the vestiges of past 
massive star formation activity. They estimated a mass 
of 1.1 x 10 5 Mq from 13 CO observations and an assump¬ 
tion of local thermodynamic equilibrium (LTE). To infer 
the total hydrogen column density, they derived a rela¬ 
tion between 13 CO and Ay, which they measured from 
star counts. Adopting a distance to G216-2.5 of 2.2 kpc 
and a diameter of 100.8 pc, Lee et al. (1994) determined 
a virial mass of 6.3 x 10 5 M 0 . They suggested that the 
factor of 6 discrepancy between the two mass estimates 
could be explained if the cloud is not in virial equilib¬ 
rium but, rather, is undergoing expansion. Nevertheless, 
even by the LTE mass estimate, G216-2.5 is more mas¬ 
sive than many other nearby Galactic, star-forming MCs 
(e.g., Imara & Blitz 2011). 

Lee et al. (1996) utilized infrared observations to in¬ 
fer the dust and low-mass star formation properties of 
G216-2.5. They derived a correlation between dust and 
gas emission reminiscent of that observed in nearby dark 
clouds. They measured a far-infrared (FIR) luminosity 
to molecular gas mass ratio of Lfir/Mh 2 < 0.7 Lq/M & . 
This ratio, which is often used as a proxy for star for¬ 
mation efficiency, is close to ~ 1 for MCs in the solar 
neighborhood (e.g., Scoville & Good 1989) and is ~ 10 
for the Rosette (Heyer et al. 2006). Lee et al. (1996) ex¬ 
amined objects in the IRAS Point Source Catalog in the 
vicinity of the cloud, and with follow-up near-infrared 
observations in the J , H , and K bands, they identified 
four of the IRAS sources as possible sites of low-mass 
star formation (Figure 1). Altogether, they estimated 
that 14 young stellar objects (YSOs) are located at these 
sites. Notably, none of these sources are located in the 
main body of the cloud. Two are associated with satellite 
clouds to the north of the main MC, and the other two 
are located near the boundary of CO emission. In this 
study, Lee et al. (1996) also pointed out that infrared 
emission in the north-east region of G216-2.5 is contami¬ 
nated by emission from the more distant H II region S286 
(Figure 1). 

A little more than a decade passed before Megeath 
et al. (2009) found evidence of low-mass star forma¬ 
tion within the main body of G216-2.5. This group took 


observations with the Spitzer Telescope’s Infrared Ar¬ 
ray Camera (IRAC) and Multiband Imaging Photome¬ 
ter for Spitzer. Using color-color selection criteria, they 
identified 74 YSOs corresponding to a peak in CO emis¬ 
sion located near (l,b) = (216.7°,—2.75°), (see Figure 
1). They noted that the surface density distribution of 
YSOs in G216-2.5 is similar to the distribution in Taurus 
and much less than the surface density of YSOs in Orion 
A. Curiously, G216-2.5 is roughly 10 times as massive as 
Taurus. 

The remainder of this work will explore the molecular 
gas distribution of G216-2.5, as revealed in rich detail by 
Herschel. I will derive and make detailed comparisons 
between the H 2 column density maps of the G216-2.5 and 
Rosette MCs, in order to probe the evolutionary state of 
the former. I will also address the question of whether 
the CO-to-H 2 conversion factor, the so-called X-factor, 
differs depending on local cloud conditions. While H 2 is 
the most abundant molecule in the interstellar medium, 
it cannot be observed directly at the cold temperatures 
of molecular clouds because it does not have a perma¬ 
nent dipole moment. Therefore, alternative tracers of 
molecular gas must be used to infer the properties of 
molecular clouds. The second most abundant molecule, 
CO, is commonly used, especially in external galaxies. 
To convert integrated CO emission, Iqo, into H 2 column 
density, X(H 2 ), one can use the X-factor, 

Xco = cm -2 (K km s -1 ) -1 . (1) 

I CO 

Bolatto et al. (2013) discussed in detail the independent 
techniques employed to determine Xco- One involves 
using tracers of column density, such as dust or opti¬ 
cally thin molecular lines. Historically, this was the first 
method used (Dickman 1975) and is justified since dust is 
observed to be well-integrated with gas (e.g., Boulanger 
et al. 1996) and since the H 2 molecule forms most ef¬ 
ficiently on dust grains. Dickman (1975) observed the 
13 CO emission from 68 Galactic clouds and compared 
this to optical extinction measurements. If a typically 
observed ratio of /( 13 CO) to /( 12 CO) is applied to those 
observations, this yields Xco = 2.2 x 10 2 °. The method 
used to derive Xco in this analysis belongs to the same 
category as that initiated by Dickman (1975). In this 
case, however, dust far-infrared emission (as opposed to 
dust extinction) is used to trace the total hydrogen distri¬ 
bution. To convert from dust to total gas to H 2 requires 
certain assumptions, including a gas-to-dust ratio (see 

§4). 

Another type of Xco determination requires the po¬ 
tentially problematic assumption that MCs are in virial 
equilibrium, and the third compares CO emission with 
7 -ray emission produced by cosmic-ray interaction with 
H 2 molecules. While there is general accord among 
the different techniques to measure Xco; variations in 
Xco have nevertheless been observed. Bolatto et al. 
(2013) propose a Galactic conversion factor of 2 x 10 2 ° 
with ± 30% uncertainty, noting that departures from 
this value should be anticipated. Dame et al. (2001) 
showed that average value of Xco hr the Solar neighbor¬ 
hood is 1.8 ± 0.3 x 10 2 ° cm -2 (K km s -1 ) -1 and varies 
as a function of Galactic latitude. The conversion fac¬ 
tor has been found to vary in spatially resolved Galactic 
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Figure 1. Integrated 12 CO emission of the G216-2.5 molecular cloud. Contour levels are from 2 to 12 K km s —1 and spaced 2 
K km s” 1 apart. The gray stars indicate the locations of IRAS point sources identified by Lee et al. (1996). The dashed circle shows the 
region where Megeath et al. (2009) identified 71 protostars and YSOs. S286 is a distant H II region unrelated to the molecular cloud. 


MCs (e.g., Pineda et al. 2008, 2010) and also in dwarf 
galaxies (e.g., Imara & Blitz 2007; Leroy et al. 2007) 
that have low metallicity with respect to the Galaxy. It 
is one of the aims of this study to address the unsettled 
questions of how and what other environmental factors 
drive changes in Xco- 

3. DATA 

To determine the quantity and distribution of molec¬ 
ular gas across the molecular clouds, I use far-infrared 
data from the Herschel Science Archives. Images in 
the PACS wavelengths (70 pm and 160 pm) and SPIRE 
wavelengths (250 pm, 350 pm, and 500 pm) that had 
passed Level 2 processing were available for regions of 
both G216-2.5 (OT2 PI: J. Kaufmann) and the Rosette 
(Motte et al. 2010). For G216-2.5, only the SPIRE im¬ 
ages are used for this analysis, since they cover the en¬ 
tire region of interest. (Given the weakness of the PACS 


emission at the low temperatures characteristic of G216- 
2.5, its exclusion is not expected to significantly modify 
the results of the fits. As will be seen in §5, the global 
properties derived for G216-2.5 are consistent with those 
arrived at in previous studies using independent meth¬ 
ods.) For the Rosette, both the PACS and SPIRE data 
are used, though they only cover a partial region (roughly 
1.5° x 1.5°) of the cloud to the east of the OB cluster NGC 
2244. 

Both clouds lie roughly two degrees below the Galactic 
midplane. G216-2.5 is centered at Galactic coordinates 
(l,b) = (216.5°, —2.5°); a distance of 2.2 kpc is adopted 
(Lee et al. 1991). The Rosette is located at ( l,b ) = 
(207°,—2°), at a distance of 1.33 kpc (Lombardi et al. 
2011 ). 

The CO data were originally obtained with the 
FCRAO 14 m telescope and were graciously provided 
by Y. Lee (originally published in Lee et al. 1994) for 
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G216-2.5, and by M. Heyer for the Rosette (see Heyer et 
al. 2006). For G216-2.5, the spatial and velocity reso¬ 
lutions of the 12 CO data cube are 50" and 0.65 km s -1 . 
For the Rosette, the spatial and velocity resolutions are 
20" and 0.127 km s -1 . The average rms noise levels in 
the 12 CO data cubes are 0.3 and 1.0 K, respectively. 

4. METHODS 

In this section, I describe the procedure used to de¬ 
termine the dust temperature, T, and column density 
distributions. 

Modified blackbody fits — The starting assumption 
is that the far-infrared emission in the MCs arises from 
cold dust which may be treated as a modified blackbody 
of the form: 

I v = B v (T)W — e~ Tv ] (2) 

where is the Planck function, given by: 

_ 2 hv 3 1 

Dv ~ c 2 e hu/kT _ l 

The dust optical depth, t u , is given by 

t v = /j,mnK v N(R) x (1/100) (4) 

where p = 1.36 is the helium correction, toh is the mass 
of a hydrogen atom, N( H) is the total hydrogen gas col¬ 
umn density, and is the dust opacity in units of cm 2 
g _1 . The factor of 1/100 in Equation 4 is the assumed 
dust-to-gas ratio. The dust opacity depends on the prop¬ 
erties and distribution of the dust and is treated as a 
power law (following Beckwith et al. 1990), 

Kv = 0,1 (lOOO GHz) ’ ^ 

where /3 ~ 2 is the dust emissivity index. 

Each pixel in a PACS or SPIRE map represents the 
flux at a certain wavelength, with three unknown values 
of interest: T, iV(H), and (or /?). The first step is 
to convolve all the images to a common resolution, that 
is, to the poorest resolution of 36", corresponding to the 
SPIRE 500 /mr data. Next, each image is projected onto 
a grid having the same number of elements, pixel size, 
and units. (The PACS and SPIRE data are originally in 
units of Jy/pixel and MJy/steradian, respectively.) To 
convert the observed flux of the images into intensity, it 
is necessary to take into account the beam size at dif¬ 
ferent wavelengths. Finally, the data are fit on a pixel- 
by-pixel basis using the MPFITFUN procedure from the 
Markwardt IDL Library (Markwardt 2009), leaving T 
and -/V(H) as free parameters and fixing the dust emis¬ 
sivity index at /3 = 2. The temperature was limited to 
the range 3 to 55 K; iV(H) was restricted to vary between 
0 and 10 24 cm -2 . 

The X-factor For each molecular cloud, the 
Xco map is determined from the ratio of the IV (H 2 ) and 
integrated 12 CO (hereafter, CO) intensity maps (Equa¬ 
tion 1). The H 2 column density is determined by sub¬ 
tracting the H I contribution from the total hydrogen 
column density, using IV(H) = IV(H I) + 2IV(H 2 ). The 
average H I column density, as measured from the LAB 
H I survey, is ~ 1.0 x 10 21 cm -2 toward G216-2.5 and 
~ 1.7 x 10 21 cm -2 toward the Rosette. The integrated 
12 CO intensity, Ico, is derived by integrating the CO 


emission over a certain velocity range. The boundaries 
of the clouds are defined by the 3 -ctrms level. For G216- 
2.5, the CO emission is integrated from 17 to 37 km s _1 , 
and the 1 -ctrms level is 0.5 K km s _1 . For Rosette, the 
velocity range is 5 to 25 km s _1 ; the I-ctrms level is 0.6 
K km s _1 . 

Since the CO and Herschel data sets have differ¬ 
ent resolutions, they are regridded onto maps having 
the same pixel size before the ratio is taken. To 
achieve similar linear resolutions at the adopted distances 
to the clouds, thereby facilitating comparison between 
them, the 7V(H 2 ) maps derived from Herschel data are 
smoothed to the resolution of the CO data 50" (= 0.53 
pc at a distance of 2.2 kpc); and the Rosette maps are 
smoothed to an angular resolution of 80" (= 0.52 pc at 
1.33 kpc). 

5. RESULTS 

5.1. Distribution of molecular gas 

Figures 2 and 3 show the H 2 column density maps de¬ 
rived from Herschel data, with CO contours overlaid. 
The average column density of H 2 gas in G216-2.5 is 
~ 7.8 x 10 20 cm -2 , and the total molecular mass is 
1.2 x 10 5 Mq, very close to the LTE mass calculated by 
Lee et al. (1994) using 13 CO measurements. The bot¬ 
tom panels of Figures 2 and 3 show that Ico closely fol¬ 
lows the distribution of H 2 gas in both molecular clouds. 
While the peaks in CO emission generally correlate with 
peaks in 1 V(H 2 ), the IV(H 2 ) maps reveal much more 
small-scale structure. Moreover, in the G216-2.5 molec¬ 
ular cloud, there are a number of small H 2 density en¬ 
hancements with little associated CO emission. 

For instance, there is a knot located at l = 215.6°, 
b = —3.4° that has a peak column density of ~ 7 x 10 21 
cm -2 , which is above the star formation “threshold” 
that has been posited by some authors (e.g., Lada et 
al. 2010; see § 6 ). However, the average density in 
this knot is ~ 1.6 x 10 21 cm -2 . At a size of roughly 
0.8 pc, this corresponds to a particle volume density of 
n ss 500 cm' 3 . There are similar high-density knots 
that have peak column densities above the star forma¬ 
tion “threshold” observed in other molecular clouds, with 
Ico ~ 2 K km s _1 located at l = 215.5°, b = —3.7° 
and l = 217.6°, b = —2.6°. By comparison, the region 
around l = 216.7°, b = —2.75°, where Megeath et al. 
(2009) identified YSOs, peaks at a column density of 
~ 2 x 10 21 cm -2 and Ico ~ 12 K km s _1 . In the satel¬ 
lite cloud located to the north-west of the main MC, 
there are knotty, filamentary-like structures that are all 
bounded by CO contours between ~ 2 and 4 K km s _1 . 
The nearly horizontal filament to the far west of the satel¬ 
lite, the location of IRAS point sources (Lee et al. 1996), 
is roughly 13 pc x 1.5 pc. It has a column density of 
~ 5 x 10 21 cm -2 and Jco ~ 10 K km s _1 . Extending 
along the eastern border of the satellite is a ~ 12 pc 
x 1.5 pc structure with 1V(H 2 ) ~ 2 x 10 21 cm -2 and 
Ico ~ 6 K km s _1 . In §5.3, these high-density regions, 
both star-forming and non-star-forming, will be exam¬ 
ined more closely. 

Within the mapped region of the Rosette (Figure 3), 
(AT(H 2 )) « 1 . 8 x 10 21 cm -2 , and the total molecular mass 
is 3.7 x 10 4 Mq. Considering that all of the molecular gas 
has not been counted, the value estimated here is in good 
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Figure 2. G216-2.5: Bottom panel - 7V(H2) map overlaid with 12 CO contours at 2, 4, 8, and 16 K km s —1 . The resolution of the maps 
are 36 /; and respectively. Top panel - Xco map. To derive Xco> the N(H. 2 ) map is convolved to the same resolution as the 12 CO 
map. 
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3. Rosette: Bottom panel: iV(H 2 ) map overlaid with 12 CO contours at 2, 4, 8, and 16 K km s 1 . The resolution of the maps 
and 80 // , respectively. Top panel - Xco map with CO contours overlaid in gray. 
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Figure 4. PDF of H 2 column densities in G216-2.5. Only pixels 
corresponding to regions where Ico > 3 — ctrms are used. The 
PDF peaks at 7V(H2) = 6.8 x 10 2() cm -2 and has a power-law tail 
whose slope in log-log space is —6.2. These and other parameters 
are summarized in Table 1. 



Log N(H 2 ) [cm" 2 ] 

Figure 5. PDF of H 2 column densities in the Rosette. Only 
pixels corresponding to regions where Iqo > 3 — ctrms are used. 
The PDF peaks at 1V(H2) = 1.4 X 10 21 cm -2 and has a power- 
law tail whose slope in log-log space is —6.3. These and other 
parameters are summarized in Table 1. 

agreement with previous determinations (e.g., Heyer et 
al. 2006; Imara & Blitz 2011). Heyer et al. (2006) found 
a mass nearly 4 times higher, but they assumed a greater 
distance to the Rosette (1600 pc) than the distance used 
here and so effectively integrated over a larger area. Us¬ 
ing this distance would increase the mass calculated here 
only by a factor of (1600/1330) 2 = 1.4. Figure 3 in¬ 
dicates that for the area mapped, there are fewer places 
(compared to G216-2.5), where there are peaks in JV(H 2 ) 
without corresponding peaks in Ico- 
To get another perspective on how high-density gas is 
distributed in the MCs, I consider the probability dis¬ 
tribution function (PDF) of the column density and the 
cumulative mass fraction (CMF) distribution. Figures 4 
and 5 display the normalized PDFs of the MCs in log-log 


Table 1: Molecular cloud PDF parameters. 


Parameter 

C216-2.5 

Rosette 

(Log iV(H 2 ))peak 

20.8 ±0.1 

24.1 ± 0.1 

A/HDpeak 

6.8 x 10 20 cm" 2 

1.4 X 10 21 cm -2 

m 

—6.2 ±0.1 

-6.3 ±0.1 

TP 

2.0 X 10 21 cm -2 

2.0 X 10 21 cm -2 


Note. (Log Ay)peak and Ay peak are the values where the 
PDF peaks, m is the slope of the high-Ay power-law tail of the 
PDF. The transition point (TP) is the value of N(H. 2 ) where 
the PDF transitions from log-normal to power-law. 


N(H 2 ) [cm- 2 ] 

10 20 10 21 10 22 



Figure 6. The normalized cumulative mass profiles as a function 
of .V f H 2 ) (equivalently, visual extinction Ay due to H 2 ) of the 
G216-2.5 and Rosette MCs. 

space, in which a log-normal distribution would appear 
as a parabola and a power-law as a straight line. In¬ 
deed, the PDFs take on the classical log-normal shape 
with a power-law tail at high column densities observed 
in other MCs (see e.g., Kainulainen et al. 2009; Froebrich 
& Rowles 2010; Lombardi et al. 2014) and predicted in 
theoretical studies (e.g., Federrath & Klessen 2013; Clark 
& Glover 2014; Ward et al. 2014). For a range of N( H 2 ) 
values before the power-law tail begins, the logarithm of 
data points are fit by a Gaussian 

| (log JV(H 2 )-fi ) 2 

/(log V(H 2 ); /i, a) = —= e ^ , (6) 

cry 27 t 

where is where the distribution peaks, effectively pro¬ 
viding the average iV(H 2 ) in the MC, and a is the dis¬ 
persion. The high-extinction portions of the PDFs are 
fit by a linear function (in log-log space). The Gaus¬ 
sian fit is robust: changing the range of data points fit¬ 
ted does not significantly alter the results, which are 
presented in Table 1. The G216-2.5 PDF peaks at 
1V(H 2 ) = 6.8 x 10 20 cm -2 and has an excess over log¬ 
normal at low column densities. The PDF of the Rosette, 
which peaks at _/V(H 2 ) = 1.4 x 10 21 cm -2 , is truncated 
at low column densities the low-density gas in the cloud 
is not fully represented in the Herschel map. Also note 
that the G216-2.5 PDF has a slightly shallower power- 
law tail than that of the Rosette. Possible implications 
of these results will be discussed in §6. 

The plot of normalized CMFs displayed in Figure 6 has 
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Table 2: Derived -Yco [10 20 cm 2 (K km s *) 1 ], 


*co 

G216-2.5 

Rosette 

Mean 

2.2 ± 1.3 

2.8 ±2.1 

ICO VS. Ay 

0.46 + AA- 
J co 

0.87 + -jlA 
laa 


horizontal axes in units of both column density and visual 
extinction, using the conversion X(H 2 )/Av = 0.94 x 10 21 
mag -1 cm -2 . The CMFs show that over a wide range of 
Ay, spanning ~ 0.5 to 15 mag, the Rosette has a higher 
fraction of its mass locked up in dense gas than does 
G216-2.5. It is noteworthy that the G216-2.5 CMF falls 
rapidly at low extinctions and then becomes shallower 
starting around 1-2 mag, revealing yet again that most 
of its mass is in the form of low-density gas. The total 
molecular mass of G216-2.5 and the mapped region of the 
Rosette are 1.2 x 10 5 Mg and 3.7 x 10 4 M 0 , respectively. 

There may be appreciable region-to-region variations 
in G216-2.5 that impact its ability to form stars or that 
are the result of past episodes of star formation. Thus, 
in §5.3, I will divide the MC into different regions and 
examine their individual PDFs. Consequences of using 
different values of /3 are discussed in 5.4. 

5.2. The X-factor 

The top panels of Figures 2 and 3 show the spatial 
variations of Xco across G216-2.5 and the Rosette. In 
both clouds, Xco is determined from Equation 1 only 
for regions where Ico is greater than the 3-orms level. 
In G216-2.5, Xco increases from ~ 0.5 x 10 20 in the cen¬ 
tral region to ~ 5 x 10 20 at the periphery. In the region 
of G216-2.5 where Megeath et al. (2009) identified em¬ 
bedded sources, Xco is slightly smaller, ~ 1.4 x 10 20 , 
than the typical Galactic value recommended by Bolatto 
et al. (2013). Figure 3 shows that for the Rosette, too, 
Xco ranges from ~ 0.5 to ~ 5 x 10 20 . Using the coor¬ 
dinates of young clusters from the Ybarra et al. (2013) 
study, YSOs in the Rosette tend to be located in regions 
that have, on average, higher values of Xco (~ 2.6x 10 20 ) 
than G216-2.5. 

The histograms in Figure 7 show that the two molecu¬ 
lar clouds have similar distributions of A"co- For G216- 
2.5, the mean value is (Xco) = (2.2 ± 1.3) x 10 20 ; for 
the Rosette, (Xco) = (2-8 ± 2.1) x 10 2 °; (see Tabic 2). 
Uncertainties in measurements of Xcoj (Tx col are de¬ 
termined by propagating the uncertainties in the X(H2) 
and CO maps. 

Another way to determine the global value of Xco is 
to take the ratio of the mean cloud properties: i.e., 
(XH 2 )/(/co)- This method is perhaps more suitable for 
galaxy-wide studies of MCs, in which the goal is often 
to determine the masses of whole clouds. For G216-2.5, 
(X CO ) = (XH 2 )/(/ CO ) = (1.8 ± 1.5) x 10 20 . For the 
Rosette, (XH 2 )/(/co) = (2.5 ± 2.4) x 10 20 . These mean 
values are somewhat lower than those cited in the above 
paragraph, since the former are more influenced by large 
values of Xco at the edges of the clouds. Neverthe¬ 
less, within the stated uncertainties, the global values of 
Xco estimated here are in agreement with one another. 

A third way to investigate Xco is to compare Ico an d 
Ay (e.g., Pineda et al. 2008; Lee et al. 2014). In Figure 
8, Ico is plotted as a function of Ay, the extinction due 
to H 2 . As noted in §4, the column density images of the 
two MCs have been regridded to maps having the same 



Figure 7. Histograms of Log(Xco)- The mean values for G216- 
2.5 and the Rosette, respectively, are 2.2 x 10 20 and 2.8 X 10 20 . 
The dotted line indicates the typical Galactic value of 2 x 10 20 . 



Figure 8. Ico versus Ay: Each data point represents values 
of a pixel in Figures 2 and 3. Overplotted is the best fit from 
Bayesian linear regression: Ay = m/co + b + e, where m is the 
slope, b is the Ay below which no CO is expected to be detected, 
and e is the intrinsic scatter about the regression. For G216-2.5, 
Ay = 0.049/co + 0-59 and e = 0.06; (equivalently, Xco — [0.46 + 
5.6//co] x 10 20 ). For the Rosette, Ay = 0.093/co + 1.0 and 
e = 0.46; (X CO = [0.87 + 9.6// C o] x 10 20 ). 

resolution as the Ico maps. The error bars in Figure 8 
are from the dispersion in each pixel of the X(H 2 ) maps. 
Error bars for Ico are not shown, but the 1-a uncertainty 
in the CO emission is taken into account in the following. 
A Bayesian linear regression analysis is performed on the 
data, according to 

Ay = to/co + b + e, (7) 

where the slope, to, is related to A'cco b is the value 
of Ay below which no CO emission is expected to be 
detected, and e is the intrinsic scatter around the re¬ 
gression. The fitted parameters are listed in Table 2. 
For G216-2.5, to = (0.049 ± 0.001) mag (K km s" 1 )" 1 
and b = (0.593 ± 0.003) mag. For the Rosette, to = 
(0.093 ± 0.001) mag (K km s" 1 )" 1 and b = (1.03 ± 0.01) 
mag. If the fit to the Rosette cloud were made over 
the same range of Ico values as for G216-2.5, (i.e., for 
Ico < 20 K km s _1 ), the resulting fit would be slightly 
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shallower. Varying the range of Ay values fitted, for ei¬ 
ther cloud, does not significantly change the slope of the 
fit. It is noteworthy that for G216-2.5, b is close to 0.5, 
the lower limit for detection in dust extinction observa¬ 
tions (e.g., Kainulainen et al. 2009), and is also near the 
column density at which atomic hydrogen begins to self 
shield from UV radiation and form H 2 . The best fit for 
b in the Rosette is somewhat high because low-column 
density regions of the cloud were not used in this analy¬ 
sis, and it most likely does not represent a true detection 
threshold. In their study of Xco in the Perseus MC, 
Lee et al. (2014) say that m x 0.94 x 10 21 is essentially 
Xco- Here, I do not exclude the functional dependence 
of Xco on Ico- Thus, for G216-2.5, the best fit parame¬ 
ters yield Xco = (0.46 + 5.6/Ico) x 10 20 . Using the aver¬ 
age CO integrated intensity across the cloud, (Ico) = 4-4 
K km s -1 , gives Xco ~ 1-7 x 10 20 . For the Rosette, 
the best fit parameters and (Ico) = 7.2 K km s -1 , yield 
X co = (0.87 + 9.6/Ico) x 10 20 « 2.2 x 10 20 . 

5.3. Regional variations in G216-2.5 

Much of the scatter in the X-factor and the Ico~ 
Ay relation is due to region-to-region variations in the 
physical environments across a given cloud. Are some 
regions more conducive to star formation than others? 
Why? Examining these parameters in individual regions 
of G216-2.5 could provide insight into the cloud’s dearth 
of star formation relative to other MCs. In particular, 
one might ask the question: do regions where there is ev¬ 
idence of low-mass star formation activity exhibit behav¬ 
ior similar to an actively, star-forming MC, such as the 
Rosette? To begin addressing this question, I select nine 
regions of G216-2.5, with the help of several plots and 
maps, for closer study. Four of the regions (1, 2, 5, and 
7) are selected due to evidence of low-mass star formation 
within them (see §2). The other regions are selected for 
having high column densities similar to the star-forming 
sites. Figure 9 re-displays the column density map of 
G216-2.5 with the regions circled and numbered. 

The PDFs of the individual regions are shown on the 
same scale in Figure 10. The histograms were generated 
from the high-resolution IV (H 2 ) map prior to smoothing; 
each plot contains all the pixels within the corresponding 
circle in Figure 9. The smallest regions have a radius of 
225" (2.4 pc) and contain an average of 812 data points. 
The size of the largest region, 1, was defined to include 
all of the embedded sources detected by Megeath et al. 
(2009); the radius is 1080" (11.5 pc), and it contains 
18,679 pixels. With the exception of region 9, each PDF 
peaks at a slightly higher column density than that of 
the PDF of the entire cloud (log X(H 2 ) = 20.8; Figure 
4), but none of them peaks at log X(H 2 ) = 21.1, where 
the Rosette PDF peaks (Figure 5). Plots 1, 2, 6, and 7 
are perhaps the best representatives of regions displaying 
some degree of log-normality with power-law tails. 

The Ico-Ay plot for each individual region is shown 
with the linear fit derived for the whole MC, for com¬ 
parison, in Figure 11. These plots are created from the 
map in Figure 2 and contain between 61 and 1440 data 
points. Regions 1, 2, and 7, where there is evidence for 
star formation, show strong evidence for 12 CO “satu¬ 
rating” at high extinctions. Interestingly regions 6 and 
3 (less so) show similar saturation, although embedded 


Tabic 3: Measurements for different values of 0. 


G216-2.5 

<N 

II 

0 = 1.8 

0 = 1.6 

(tV(Ho)) [10 2U cm' 2 ] 

7.8 

5.0 

2.8 

(Tdust) [K] 

13.2 

14.1 

15.1 

(X co > [10 2 °] 

2.2 

1.4 

0.7 

ATmc [10 5 Mq\ 

1.2 

0.78 

0.43 

Rosette 

CN 

II 

sa 

II 

bo 

CD 

II 

(N( Ho)) [10 2U cm - "] 

18 

12 

7.1 

(Tdust) [K] 

14.2 

15.3 

16.6 

(X co ) [10 2 °] 

2.8 

1.8 

1.0 

ATmc [10 5 AT©] 

0.36 

0.24 

0.15 


sources have not yet been observed to be forming in the 
dense gas there (see §6 for more discussion). Another 
noteworthy feature is the near-verticality of the plots cor¬ 
responding to regions 4, 5, and 8. In terms of Equation 7, 
this corresponds to a slope equal to zero; in other words, 
these non-star-forming regions appear to have a constant 
X-factor. 

5.4. Varying (3 and uncertainties 

The major source of uncertainty in this analysis arises 
from the modified blackbody fitting described in §4. Sin¬ 
gle values of X(H 2 ), T , or the emissivity index (3 along 
each line-of-sight do not reflect reality. Apparent vari¬ 
ations in temperature may actually be variations in /3. 
Using different values of (3 in the the modified black- 
body fitting would result in different outcomes of the 
X(H 2 ) maps and, thus, the measurements of Xco- As 
f3 decreases, dust temperature increases and X(H 2 ) de¬ 
creases, resulting in a smaller overall X-factor and less 
massive clouds. 

In this analysis, after testing different values of (3, 1 use 
a fixed value of f3 = 2 and found this to be a valid choice, 
since it yields average MC column densities and total 
masses in agreement with independent methods of mea¬ 
suring these quantities. For instance, if I had adopted a 
value of (3 = 1.8, the average temperature in each MC 
would be ~ 1 K higher; the typical column density and 
total mass would decrease by a factor of ~ 1.5. Table 
3 lists total cloud masses and typical values of X(H 2 ), 
dust temperature, and Xco measured for three different 
values of (3. 

For the G216-2.5 cloud, another source of uncertainty 
arises from the non-inclusion of the PACS data. But 
given the low dust temperatures characteristic of G216- 
2.5 ((T dust ) = 13.2 K; Table 3) and the weakness of the 
PACS flux, it is not expected that the addition of the 
PACS data would significantly modify the results. The 
dust temperature histograms of the two clouds are dis¬ 
played in Figure 12. 

6. DISCUSSION 

This analysis has measured Xco in two molecular 
clouds at sub-parsec resolutions. The mean values are 
consistent with the typical Xco ~ 2 x 10 20 measured 
in other Milky Way MCs (e.g., Bolatto et al. 2013). 
In Perseus, a molecular cloud having a similar mass to 
and perhaps more YSOs than the Rosette, Pineda et 
al. (2008) measured variations in the X-factor, Xco ~ 
(0.9 — 1.8) x 10 20 , between six subregions of the cloud. 
These authors also found that the measure of Xco is 
dependent upon saturation of the CO emission above 
Ay > 4 mag. In another study of Perseus, Lee et al. 
(2014) measured an Xco of ~ 3 x 10 19 , an order of mag- 




218 217 216 215 214 

Galactic longitude [degrees] 

Figure 9. High-column-density regions selected for closer examination. The resolution of the 7V(H2) map is 36”, the same as Figure 
2. The dashed circles circumscribe regions where star formation has been identified. Solid circles show high-.Yfllo) regions where star 
formation has not been detected. All of the pixels within each circle are included in the corresponding plots shown in Figures 10 and 11. 


nitude smaller than the mean values measured for the 
Rosette and G216-2.5 in this study. Lee et al. (2014) 
also found variations across the cloud. In their analy¬ 
sis of the Taurus Molecular Cloud, Pineda et al. (2010) 
derived Xqo ~ 2.1 x 10 20 with variations in the range 
(1.6 — 12) x 10 20 . Intriguingly, at least judging from 
the few studies of resolved MCs, it does not appear that 
Xco depends on a cloud’s evolutionary status. 

More insight into the evolutionary state of G216- 
2.5 comes from consideration of its estimated star for¬ 
mation rate and its dense gas distribution. Evidence for 
large variations in the star formation activity of molecu¬ 
lar clouds has been known for some time (e.g., Mooney 
& Solomon 1998), and G216-2.5 is no longer the sole 
example of a massive MC with a low SFR. In their se¬ 
ries of extinction mapping surveys, Lada et al. (2009) 
identified the California Molecular Cloud, a nearby (0.5 
kpc), massive (~ 10 5 M 0 ) molecular cloud with an SFR 
of ~ 70 Mq Myr -1 . In a study of 11 molecular clouds 
from their survey, Lada et al. (2010) argued that the 


SFR is determined by the total amount of gas above a 
threshold density. In terms of extinction, they estimated 
this threshold to be Ak = 0.8 ± 0.2 mag, corresponding 
to Ay = 7.3 ±1.8 mag. They compared the number of 
YSOs, Nyso, hi the clouds to their masses above the ex¬ 
tinction threshold, Mo, and found a nearly linear relation 
of the form Nyso = Mq, with a = 0.96 ± 0.11. From 
the column density maps, the total molecular mass in 
G216-2.5 having extinctions above 7.3 mag is 874 Mq, 
about 400 Mq less than what Lada et al. (2010) esti¬ 
mated for the Ophiuchus Molecular Cloud (1766 M 0 ). 
Including all of the data available on G216-2.5, the cloud 
has roughly 88 YSOs, counting the 14 identified by Lee 
et al. (1996) and the 74 of Megeath et al. (2009). For 
the Rosette, M 0 = 1389 Mq, perhaps a slight underesti¬ 
mate of the amount of dense gas, since the entire cloud 
was not mapped in this study. The known number of 
YSOs in this MC is ~ 461 (Ybarra et al. 2013). These 
quantities are plotted in Figure 13, along with the data 
from the Lada et al. (2010) sample. 
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Figure 10. PDFs of regions indicated in Figure 9. Regions where star formation activity has been identified are labeled (SF). For reference, 
the dotted line indicates the location of the peak of the PDF of the entire cloud, (LogA^(H2)) pea k — 20.8. 


For the clouds in their sample, Lada et al. (2010) as¬ 
sume the median mass for the stellar initial mass function 
(IMF) to be 0.5 M 0 (e.g., Muench et al. 2007), and they 
derive an expression relating the SFR to Nyso- 

SFR = 0.251Vyso M q Myr -1 . (8) 

Using this formulation, the SFR in the Rosette is 115M 0 
Myr -1 , roughly halfway between the rates in Taurus 
(84M 0 Myr -1 ) and Perseus (15OM 0 Myr -1 ). In G216- 
2.5, SFR « 22 Mq Myr -1 , similar to the rates in the 
Lupus 3 and RCrA molecular clouds. This corresponds 
to 44 stars Myr -1 , not far from the estimate of Megeath 
et al. (2009), who determined that G216-2.5 is forming 
ss 60 stars Myr -1 . 


Yet considering that Equation 8 assumes a star forma¬ 
tion duration of 2 Myr and that G216-2.5 has a partic¬ 
ularly low ratio of Class II to Class I YSOs, compared 
to nearby Galactic molecular clouds (Gutermuth et al. 
2009; 2011), it is possible that a higher SFR in G216- 
2.5 may be more true to reality. Given the latest ob¬ 
servations of Megeath et al. (2009) that G216-2.5 has 
33 Class I sources, which may have been born as re¬ 
cently as ~ 0.5 Myr ago, (e.g., Hatchell et al. 2007; Bat- 
tersby et al. 2014), Equation 8 could be reformulated as: 
SFR = 0.25Vp roto Mq Myr -1 (2/f proto ), where ^proto ^ 
the typical protostar lifetime in units of Myr, and Np Ioto 
is the number of protostars. In this case, the SFR in 
G216-2.5 would be ~ 33 M 0 Myr -1 . 
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Figure 11. Iqq vs. Ay plots of regions indicated in Figure 9. Regions where star formation activity has been identified are labeled (SF). 
Overplotted for reference is the best fit to the data for the entire cloud (black line). 


Lada et al. (2010) also give a predicted expression 
for the SFR that scales with dense mass: SFR « 4.6 x 
10 -2 Mo M 0 Myr -1 , where Mo is the mass of gas hav¬ 
ing Ay > 7.3 ±1.8 mag. As justification, they ar¬ 
gue that high-extinction material corresponds to material 
with high volume number density (n(H 2 ) > 10 4 cm -3 ). 
Using this equation, the predicted SFR in G216-2.5 is 
~ 40 Mq Myr -1 , about half the SFRs in Ophiuchus and 
Taurus, and consistent with the above estimates from 
Equation 8. 

If anything, the predicted value may slightly overes¬ 
timate the SFR in G216-2.5, however, due to the as¬ 
sumption that high-extinction material corresponds to 


n(H 2 ) > 10 4 cm -3 . As Clark & Glover (2014) pointed 
out, the appearance of high-density gas at high extinc¬ 
tions does not necessarily mean high-extinction material 
is composed only of high-density gas. In their simula¬ 
tions of clouds of different masses, Clark & Glover (2014) 
found that early on in a cloud’s evolution, the quantity of 
gas above the volume density threshold can be less than 
the quantity of gas above the column density threshold 
by orders of magnitude. Indeed, this may be the case in 
G216-2.5 if it is a young MC. As was shown in §5.1, G216- 
2.5 has a number of high-column-density regions with low 
number densities and with no evidence of star formation 
activity. This is in keeping with the Williams & Blitz 


















13 



10 12 14 16 18 

T [K] 


Figure 12. Histograms of dust temperatures derived from Her- 
schel data. The G216-2.5 histogram (black) peaks at 13.2 K; the 
Rosette histogram (red) peaks at 14.2 K. 



Figure 13. Plot of the number of YSOs, ATysOj vs - MC mass 
above the extinction Ak = 0.8 mag (Ay = 7.3 mag), adapted 
from Lada et al. (2010). The circles are data points taken from 
Lada et al. (2010), as is the best fit to the data (dashed line). 
The stars represent data points for the Rosette and G216-2.5. The 
filled circle represents the California Molecular Cloud. 

(1998) study of clumps in G216-2.5 and the Rosette in 
the CO(3-2), CO(l-O), and CS(2-1) lines. They found 
that ~ 20% of the clumps in the Rosette are star-forming 
or incipient star-forming and ~ 80% are “dormant” (i.e., 
non-star-forming). In G216-2.5, by comparison, nearly 
all the clumps are dormant. Moreover, Williams & Blitz 
(1998) found that dormant clumps in the Rosette tend 
to be denser and have steeper radial density profiles, sug¬ 
gesting that they are more evolved than dormant clumps 
in G216-2.5. Thus, as G216-2.5 as a whole ages, per¬ 
haps the number of gravitationally bound cores capable 
of forming stars will increase, as was the case for the 
clouds in the Clark & Glover (2014) simulations. With 
regards to the star formation history, with the data avail¬ 
able here and the given assumptions, the estimates of the 


SFR in G216-2.5 are within only a factor of 2 of each 
other. New observations in other regions of the cloud 
of deeply embedded sources may help to provide a more 
accurate measure of its star formation history. 

The idea that different regions within G216-2.5 may 
be evolving at different rates may also explain, in part, 
the regional differences in PDFs. From Figure 10, the 
following points are of note. First, the PDFs of star¬ 
forming regions display a power-law tail. But at least 
three regions in which star formation has not been de¬ 
tected (6, 8, and 9) also have PDFs with an excess of high 
column densities (relative to log-normal), leading one to 
speculate that these regions might have deeply embedded 
sources or that star formation is incipient. Second, there 
does not appear to be a characteristic width of PDFs of 
star-forming regions. While signatures of star formation 
activity have been identified in Regions 1, 2, 5, and 7, 
the widths of PDFs of Regions 1 and 7 are broader than 
those of 2 and 5. Moreover, the width of these latter two 
regions are also narrower than 3, 6, and 9, in which star 
formation has yet to be detected. Third, each of the star¬ 
forming regions have PDFs which peak at higher column 
densities than the PDF of the cloud as a whole. Fourth, 
the PDFs of regions without detected star formation 
Regions 3, 4, and 6—also peak at high column densities, 
relative to the peak of the entire cloud (indicated by a 
dotted line in Figure 10). Fifth, in addition to possess¬ 
ing high-column density gas, regions with star formation 
contain a significant amount of low-colunm-density gas, 
to varying degrees depending on the region. The different 
sizes of the regions analyzed (Figure 9) do not account 
for these regional differences, since the characteristics of 
the PDFs listed above are robust to variations in the 
boundary sizes. Factors, such as turbulence, that may 
contribute to the regional differences will be explored in 
a future study. 

Theoretical and observational studies suggest that, as 
a consequence of the gravitational collapse occurring on 
small scales within molecular clouds, the column density 
PDFs evolve. Numerical simulations of molecular cloud 
evolution indicate that the PDF of a cloud acquires a 
power-law tail or “wing” at high extinctions and that 
this tail tends to flatten over time as the star forma¬ 
tion efficiency increases (e.g., Federrath & Klessen 2013; 
Clark & Glover 2014; Ward et al. 2014). In their study of 
23 Galactic molecular clouds, Kainulainen et al. (2009) 
found that actively star-forming clouds have PDFs with 
non-log-normal tails at high-extinctions, while quiescent, 
non-star-forming clouds have log-normal PDFs over the 
entire range of column densities. Figures 4 and 5 show 
that both G216-2.5 and the Rosette have high-extinction 
tails. Interestingly, the slope of the G216-2.5 PDF is 
slightly flatter than that of the Rosette, even when the 
range of values chosen for the fits are changed. Tak¬ 
ing this to suggest that G216-2.5 is more evolved than 
the Rosette would be an over-interpretation, however, 
since it should be kept in mind that the slope of the 
high-extinction tail is probably not universal for clouds 
of similar age (Federrath & Klessen 2013). Moreover, 
observational column density PDFs do not easily trans¬ 
late to volumetric PDFs, and as discussed above, high 
AT(H 2 ) does not necessarily imply high n(H 2 ). To better 
understand (i) the regional differences in PDFs discussed 
above, and (2) how the properties of gas in G216-2.5 de- 
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termine its SFR, high-resolution observations of molec¬ 
ular lines that trace high-density gas would be useful. 
With observations of NH 3 , N 2 H, and HCN, for instance, 
future work could focus on comparing cores in G216- 
2.5 with those in other MCs. 

In the context of observational studies that aim to un¬ 
derstand the nature of diffuse gas in molecular clouds, 
the large excess over log-normal of low-colunm-density 
gas in G216-2.5 is interesting. Previous studies that de¬ 
rived column density PDFs of MCs from extinction mea¬ 
surements are sometimes inconclusive about the origin 
of this low-extinction material (e.g., Kainulainen et al. 
2009; Lombardi et al. 2011). In particular, it can be 
difficult to ascertain whether low-extinction material is 
real (i.e., from the cloud of interest) or if it is mostly 
contamination from foreground objects along the line of 
sight toward the cloud. The work presented here indi¬ 
cates that this low-extinction feature of the PDF persists 
when column densities are measured using dust emission 
measurements and, in the case of G216-2.5, is likely to 
represent real structure. 

Spatially, the low-density tail of the G216-2.5 PDF cor¬ 
responds to low-density gas from within the cloud itself. 
In the PDF represented in Figure 4, each column density 
is the excess above the boundary of the cloud, as defined 
by Ico (§3). Thus, the low-column density tail of the 
PDF does not arise from the adjacent photodissociation 
region, discussed in Williams & Maddalena (1996), which 
lies beyond the CO boundary of the cloud. Not including 
gas beyond this boundary also ensures that high-density 
features are not over-represented in the tail of the PDF. 

It is worth considering whether this low-density gas 
truly arises from within G216-2.5 itself, or is due to line of 
sight (LOS) contamination. One way to test for the lat¬ 
ter possibility is to consider how much of a contribution 
the low-density gas makes to the total mass of G216-2.5. 
For instance, gas having values to the left of the peak in 
the PDF in Figure 4 (i.e., AT(H2) < 6.8 x 10 20 cm -2 ) con¬ 
tributes to 30% of the estimated total molecular mass. 
Subtracting this contribution would result in an unreal¬ 
istically low value for the total mass, so it is likely that 
this low-density material results from real dust emis¬ 
sion within G216-2.5. Schneider et al. (2014), unlike 
this study, make a correction for LOS contamination to¬ 
ward G216-2.5, which brings down their average column 
density before correction; they consider only extinctions 
above Ay > 1 mag. In this study, the contribution of 
atomic gas to the total hydrogen column density is cor¬ 
rected for (§4), which also has the effect of driving down 
the value derived here for the average column density, 
and values below Ay = 1 mag are included in the analy¬ 
sis. Ultimately, the combined effects of our different anal¬ 
yses result in Schneider et al. (2014), deriving a higher 
value for the peak of the G216-2.5 PDF ( Ay = 1.9 mag, 
after correction) than that derived here. 

It is also instructive to consider the column density at 
which the PDF transitions from log-normal to the high- 
density power-law tail (called TP in Table 1). In terms 
of visual extinction, this transition occurs at ~ 2 mag, in 
both clouds. Kainulainen et al. (2011) called this tran¬ 
sition extinction Ay 1 ' and found that it occurs between 
2-4 mag in most of the clouds from their 2009 study, 
in which they derived column densities from extinction 


maps using star counts. Schneider et al. (2014) derived 
column density PDFs for a sample of clouds, including 
G216-2.5, from Herschel dust emission maps. For G216- 
2.5, they found that Ay 11 ss 4 mag; and overall, for their 
group of four MCs, they found a higher range of transi¬ 
tion extinctions occurring between 4-5 mag. 

In their study of cluster formation in the Rosette, 
Schneider et al. (2012) derived PDFs in six subregions of 
the cloud. They found that the PDF for the whole cloud 
has a log-normal distribution up to a visual extinction of 
~ 3 mag and a power-law tail at high extinctions (Ay = 3 
- 20 mag). In this study, the contribution of H I to the to¬ 
tal gas column density is subtracted, thus, lower values of 
the PDF peak and the transition extinction are derived: 
Ay 11 « 2.1 mag (equivalently, 1 V(H 2 ) « 2.0 x 10 20 cm -2 ; 
see Table 1). The essential point, nevertheless, is that 
the shape of the PDFs derived here—log-normal at low 
extinctions with power-law tails at higher extinctions—is 
consistent with previous studies, in spite of the different 
criteria for defining clouds. 

Taken together, the findings of this and previous stud¬ 
ies give the impression that G216-2.5 is a “young” MC 
no more than a few million-years-old. Hartmann et al. 
( 2001 ) noted that the stellar populations of most molec¬ 
ular clouds have age spreads of no more than 1-3 Myr. 
This fact as well as the absence of older stars favor MCs 
being transient entities that live no more than ~ 30 Myr 
(e.g., Blitz & Shu 1980). Hartmann et al. (2001) argue 
for a picture in which star formation and MC formation 
are triggered by the same large-scale flow of distant star¬ 
forming sites. Outflows from high-mass stars ultimately 
disperse the cloud of which they formed; but it is an 
open question as to whether low-mass stars can similarly 
disperse their parent cloud on short timescales. 

This line of reasoning is not necessarily at odds with 
the argument of Lee et al. (1994) that G216-2.5 is a rem¬ 
nant cloud from past massive star formation. One need 
only (re)consider the term “remnant” and its connota¬ 
tion of oldness. If MCs form at the interface of large- 
scale flows, as some simulations suggest (e.g., Hartmann 
et al. 2001; Hennebelle et al. 2008), then new clouds 
naturally condense from the “remnants” of old clouds. 
Future work related to this study will examine the dy¬ 
namics of G216-2.5 and of the atomic gas associated with 
it. Williams & Maddalena (1996) already showed that 
G216-2.5 appears to be connected to S287, a distant star¬ 
forming region, by an H I cloud. And Lee et al. (1994) 
established that G216-2.5 has a global velocity gradient 
that could be explained if the MC is part of a large, ex¬ 
panding shell. Another interpretation, however, is that 
the gradient is due to shear (if, for instance, the cloud 
formed at the interface of converging flows). Considered 
together, this could mean G216-2.5 recently merged from 
large-scale flows and has not had enough time to build 
up the H I necessary to shield H 2 and CO from disso¬ 
ciating UV radiation as effectively as older MCs. This 
is not a firm conclusion, nevertheless, and so subsequent 
work will concentrate on understanding the large scale 
surroundings of G216-2.5. 
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7. CONCLUDING REMARKS 

Although evidence for differing SFRs in MCs is not 
new, G216-2.5 still seems to be exceptional, considering 
the peculiar combination of its properties: (1) It has a 
total mass comparable to the most massive MCs (e.g., 
California, Orion B); (2) its total dense mass is compa¬ 
rable to medium-sized MCs with hundreds of YSOs (e.g., 
Taurus, Ophiuchus); and yet (3) it forms stars at a rate 
similar to low-mass clouds (e.g., Lupus). In this work, 
I compare the properties and distribution of molecular 
gas of G216-2.5 and the Rosette, two MCs having ex¬ 
tremely different SFRs. I use archival Herschel data to 
create high-resolution 7V(H 2 ) maps and compare these 
with 12 CO data to derive images of Xco at ~ 0.5 pc 
spatial resolution for both clouds. The main results are 
as follows: 

1 . For G216-2.5, (X C o) = 2.2 x 10 20 , and for the 
Rosette, (Xco) = 2.8 x 10 20 , values consistent with 
the typical Aco measured in other Milky Way MCs 
(e.g., Dame et al. 2001; Bolatto et al. 2013). The 
conversion factor increases to as high as ~ 5 x 10 20 
toward the edges of the MCs but hovers between 
~ 0.5 and 2 in the main bodies of the clouds. 

2. The AT(H 2 ) PDFs of the clouds both have log¬ 
normal distributions with power-law tails at high 
column densities. The G216-2.5 PDF peaks at 
N( H 2 ) = 6.8 x 10 20 cm -2 (Ay = 0.6 mag); that 
of the Rosette peaks at 1V(H 2 ) = 1.4 x 10 21 cm -2 
(Ay = 1.4 mag). The Rosette has a higher frac¬ 
tion of its mass locked up in gas having extinctions 
ranging from Ay ~ 0.5 to 15 mag. 

3. To compare the properties of sites with star for¬ 
mation and high-density areas without stars, nine 
regions in G216-2.5 are are selected for closer ex¬ 
amination. The PDFs of the star-forming regions 
and at least two of regions where star formation 
has not been observed appear to display power- 
law tails, a feature that previous authors have at¬ 
tributed to the dense gas associated with star for¬ 
mation. Eight of the regions peak at higher extinc¬ 
tions than the cloud as a whole, but there is not a 
significant difference in the peak location between 
star- and non-star-forming regions. Also, none of 
the regions peak close to Ay peak measured in the 
Rosette. 

4. According to the Ico-Ay plots of the different 
regions, there is evidence for CO saturation at 
high extinctions in both star-forming and non-star- 
forming regions. Additionally, while in some re¬ 
gions Aco appears to be a function of Icq > in oth¬ 
ers, Aco is nearly constant. 

5. The Rosette has a higher fraction of its mass in 
the form of dense gas and contains 1389M 0 of gas 
above the so-called extinction threshold for star for¬ 
mation, Ay = 7.3 mag. The G216-2.5 cloud has 
874M 0 of dense gas above this threshold. 
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